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Abstract 

We consider two independent bosonic oscillators immersed in a common 
bath, evolving in time with a completely positive, markovian, quasi- 
free (Gaussian) reduced dynamics. We show that an initially separated 
Gaussian state can become entangled as a result of a purely noisy mech- 
anism. In certain cases, the dissipative dynamics allows the persistence 
of these bath induced quantum correlations even in the asymptotic equi- 
librium state. 
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1. INTRODUCTION 



When two non-interacting, independent systems are immersed in a common bath, 
decoherence effects are expected to arise counteracting any quantum correlation initially 
present among the two subsystems [1-5]. In certain circumstances however, the environ- 
ment can also enhance entanglement. Most commonly, this phenomenon is the result of 
an hamiltonian coupling between the two subsystems [6-10] generated by the action of the 
bath; remarkably, it can also occur in the markovian regime through a non-hamiltonian, 
purely noisy mechanism [11-22]. 

This environment induced entanglement generation has been identified so far in models 
involving finite-level systems. In this note we shall instead examine the case of simple 
infinite dimensional systems following a purely markovian dynamics. More specifically, we 
shall study the behaviour of two independent bosonic oscillators evolving in time according 
to a semigroup of completely positive maps. We shall limit our considerations to the 
phenomenologically relevant set of quasi-free (or Gaussian) states [23-29] and to dynamics 
that preserve this set, the so-called quasi- free semigroups, whose properties are well studied 
in the literature [30, 31]. 

After a brief review of the theory of quasi-free dynamical semigroups in relation to 
the present situation, we shall first discuss the possibility of entanglement creation at the 
beginning of the evolution, in the vicinity of the initial time t = 0. Using the partial 
transposition criterion for two-mode Gaussian states, we shall derive the condition that 
assures environment induced entanglement generation as soon as t > 0. 

Although in general entanglement will keep growing also away from the neighborhood 
oft = 0, the decoherent, noisy character of the environment will eventually set in. Whether 
or not quantum correlations are actually left in the asymptotic long time regime is the result 
of a delicate balance between entanglement generation and decoherence. It is remarkable 
that for certain environments the long time equilibrium state could indeed retain a non- 
vanishing entanglement, and an example of this phenomenon is explicitly presented at 
the end. In view of the direct phenomenological relevance of Gaussian states in quantum 
optics, we believe that these results might be of help in the design and actual realization 
of quantum devices and circuits. 



2. GAUSSIAN STATES 

As explained in the introductory remarks, we shall study the dynamics of independent 
oscillators in weak interaction with a large environment, so that their reduced time evo- 
lution can be well represented by a markovian, completely positive dynamics. Although 
our considerations concerning bath-induced entanglement creation will involve two such 
oscillators, for sake of generality in this Section and the next will let their number n to be 
arbitrary. The states of this system will be represented by a density matrix p, i.e. by a 
positive hermitian operator, with unit trace, acting on the bosonic Hilbert space Ti. It can 
be identified with the tensor product of n independent Fock spaces; they are generated 
from the vacuum state through the action of polynomials in the creation a\ and annihila- 
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tion aj, i = 1, 2, . . . , n, operators pertaining to the single oscillators. These operators obey 
the standard bosonic oscillator algebra: [c^, aj] = 8^, [a}, aj] = = [aj, a,]. 

In dealing with this algebra, it is useful to adopt the holomorphic representation 
[32-35]; it allows working with explicit expressions for both states and operators, including 
the density matrix p. In this formulation, the elements of the bosonic Hilbert space 
H are represented by holomorphic functions ip(z) of the set of n complex variables z = 
(z\,Z2, ■ ■ ■ ,z n ), with inner product :^ 

= f r(z)Hz) e"^^ Uidzidzi , (2.1) 

where * signifies complex conjugation. To every operator O acting on 7i there corre- 
sponds a kernel 0(z,z) of 2n independent complex variables z = (zi, z%, . . . , z n ) and 
z = (z\, Z2-, . . . , z n ), such that for the state \<f>) = 0\ip) one finds the representation: 

0(z, w) ip(w) e~ 51i mm n, dw z dw z . (2.2) 

In particular, the creation and annihilation operators, when acting on a state \ifj), are 
realized by multiplication and differentiation by the variables zf. 

a\\iP) -> z^(z) , a#) -> ^(z) , (2.3) 

while the identity operator is represented by exp(^ ZiZi). 

Among all density matrices p for the n oscillators, the so-called quasi-free or Gaussian 
states are of particular interest: they can be easily produced in experiments in quan- 
tum optics [27-29]. In the holomorphic formulation, they can be defined as possessing a 
representation kernel p(z, z) of generic Gaussian form: 



p(z, z) = VAf exp 



- Z T • G 1 ■ Z + Z i Z i 
i 



(2.4) 



where z T is a 2n-dimensional row vector with components (z\, Z2, ■ ■ ■ , z n , zi, Z2, ■ ■ ■ , z n ), 
while z is its corresponding column vector (T signifies matrix transposition and • matrix 
multiplication). The covariance G is 2n x 2n matrix; with respect to the natural decom- 
position z T = (z z) into the set of variables Zi and Zi, it can be conveniently parametrized 

as 

G=(t I) , (2-5) 



1 a* 

where a and (3 are n x n matrices, with a symmetric and (3 hermitian to assure the 
hermiticity of p. They represent the averages in the given state p of quadratic operators 



t Here and in the following we use the conventions of Ref. [32] 
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in the creation and annihilation operators: 

(a< aj) = Tr [en aj p]= J q- q- [p(z, z)] e ^ ZkZk U k dz k dz k = ay , 

(aj at) = Tr [aj a] p] = J z z d p(z, z) e" ^ ~ ZkZk U k dz k dz k = a* tJ , (2.6) 

/d \ ^ 

— [% p{z, z)] e" 2feZfe n fc d^fc d^ = /% • 

As a result of the last relation, the matrix f3 turns out to be non-negative. For simplicity, 
in writing (2.4) we have assumed (aj) = (aj) = 0; this condition can be easily released 
starting with a more general Ansatz for p(z,z), containing in the exponential also linear 
terms in Zi and zf, it will not be needed for the considerations that follow. The trace 
condition, 

Tr[p] = J p(z, z) e~ ~ ZiZi U t dz t d Zi = 1 , (2.7) 

further fixes the normalization constant, M = det(G), provided the previous integral 
makes sense [32]. Indeed, convergence of (2.7), as well as of the integrals in (2.6), put 
further constraints on the entries of the covariance matrix G, or equivalently on those 
of a and (3. As similarly done with the vector z, let us collect the 2n annihilation and 
creation operators aj and aj into the column vector a and its hermitan conjugate row vector 
ee (a| aj) = (a|, a|, . . . , a^, ai, a2, . . . , a n ). Then, necessarily, the following 2n x 2n 
matrix of bilinear expectation values results non negative: 

(a M at)EETr[a M at p] >0, p, v = 1, 2, . . . , 2n . (2.8) 

It turns out that (2.8) is also a sufficient condition for the expression in (2.4) to represent 
a physical state [23, 36-38]. 

It is customary to rewrite this condition in terms of the 2n x 2n matrix V of symmetric 
bilinears: 

V^ = ^({a M , 4})^Tr[( aAt at +ata M )p] ; (2.9) 



2X i ^ V j/ 2 

in terms of the natural decomposition of a^ = (aj aj) into the set of all creation and 
annihilation operators, it can be explicitly written as: 

V=(£ It)-\- (2-10) 



d* 1 J 2 
Introducing also the 2n x 2n matrix of commutators 



[ a M, 4])= (j , (2.11) 



one can finally express the positivity condition (2.8) as: 



V + — >0. (2.12) 



In the following, we shall limit our considerations to the set of Gaussian states, i.e. to 
the states p represented by kernels of the form (2.4) satisfying the condition (2.8), or 
equivalent ly (2.12). 

A particularly important class of Gaussian states are the pure ones. In the holomor- 
phic representation they are described by properly normalized Gaussian functions of the 
variables zf. 



y n (z) = det 1/4 (l - 0*0) exp -- ^ z t Q 



(2.13) 



with O a complex, symmetric matrix, such that |0| < 1 to guarantee norm convergence. 
The corresponding kernel p{z,z) = ^n(z) ^q(z) can be cast in the form (2.4), with sub- 
matrix coefficients a = —0(1 — 0*0) _1 and $ = (1 — 0*0) _1 , respectively. 



3. QUASI-FREE QUANTUM DYNAMICAL SEMIGROUPS 

Our analysis is based on the assumption that the time evolution of the set of inde- 
pendent oscillators immersed in the common bath be markovian and given by a quantum 
dynamical semigroup; this is a completely positive, trace preserving, one parameter family 
of linear maps, acting on the set of density matrices p representing the oscillator states. 
These maps are generated by equations of the following Kossakowski-Lindblad form [39-41] : 

^ =C[p(t)] ee -i[H,p(t)]+L[p(t)] , (3.1) 

with H and effective hamiltonian and L a dissipative piece, that can be abstractly written 
as 

L[p] = Wl* p4 - \{L\L ki p]\ (3.2) 

k ^ ~ ' 

in terms of a collection of suitable, well-behaved operators Lf~. 

Being interested in the set of Gaussian states, we would like to characterize those 
quantum dynamical semigroups that preserve that set. Although already discussed using 
abstract, mathematically rigorous techniques [30, 31], this problem has a simple and direct 
solution in the holomorphic representation: the semigroup generated by (3.1) leaves the 
form (2.4) for p invariant provided Lk is a linear and H a quadratic combination of the 
creation and annihilation operators. Indeed, only in this case the r.h.s. of (3.1) results 
quadratic in the variables Zi and their derivatives when acting on the representation kernel 
p(z,z), thus preserving its generic form (2.4); as a consequence, the evolution equation 
(3.1) reduces to a linear, differential equation for the entries a and j3 of the covariance 
matrix G (given explicitly in (3.7) below). 

In view of this result, the effective hamiltonian will be taken to have the generic 
quadratic form 

1 n 
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whith the matrix Co hermitian and positive. Note that terms containing a^aj and their 
hermitian conjugate can be eliminated by a suitable unitary canonical transformation 
[32, 34, 35], and therefore do not appear in (3.2). Similarly, the most general allowed 
dissipative term takes the following form: 



1 n f 

L [p] = g J2\ *li3 ( [ a iP> a l] + l a v P a l] ) + &ij ( [ a ]^' ai \ + [°]' ^] ) 



(3.4) 



where the n x n coefficients matrices ry, a and A encode the physical properties of the 
environment and can be expressed in terms of the Fourier transform of the correlation 
functions in the bath [1-5]. 

These parameters are not completely arbitrary. First of all, fj and & need to be 
hermitian in order to comply with the hermiticity preserving requirement of the generated 
semigroup. In addition, the request of complete positivity gives rise to further constraints. 
Indeed, by using the 2n-dimensional vectors a and at introduced in the previous section, 
the dissipative term in (3.4) can be recast in compact form as: 

L \p\ = Yl c ^( a vp a i- 2-{ a i a ^ p}) ' ( 3 - 5 ) 

where the bath coefficients are now embedded in the 2n x In Kossakowski matrix C. It 
is well known that the requirement of complete positivity of the dynamics generated by a 
dissipative term L[p] in the form (3.5) is equivalent to the positivity of the Kossakowski 
matrix [39-41]; in the present case, using the block decomposition introduced before in 
writing (2.10), this condition explicitly reads: 

C=(f^)>0. ,3.6) 

The finite-time evolution maps obtained from the equation (3.1), with hamiltonian as 
in (3.3) and dissipative term as in (3.5), with C > 0, are known in the literature as quasi- 
free quantum dynamical semigroups [30, 31]: as already observed, they are characterized 
by the property of transforming the set of quasi- free (Gaussian) states into itself. 

Inserting the general Ansatz (2.4) for the kernel p(z, z; t) in the evolution equation 
(3.1), with (3.3) and (3.4), and using the prescriptions (2.3), one easily derives the equation 
obeyed by the time-dependent covariance G(t), or equivalently those for the sub matrices a 
and (3. It can be more conveniently rewritten as an equation for the symmetric covariance 
matrix V introduced in (2.9): 

dtV(t) = A* • V(£) + V(t) ■ A + B . (3.7) 
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The 2n x 2n matrices A and B contain the dependence on the hamiltonian i2>ij and dissi- 
pative coefficients rjij, d~ij, A^-; using the same decomposition introduced in the definition 
of V in (2.9), one explicitly finds: 






-2(A( S ))* 



a + fj* 



(3.8) 



~ (S) " (A) 

where A^- and A^- are respectively the symmetric and antisymmetric parts of the complex 
matrix \j. 

The evolution equation (3.7) has been originally derived using the equivalent Heisen- 
berg picture where the time-evolution affects the observables O instead of the states p; the 
two pictures are connected by the duality relation involving expectation values 



(0)(t) = Tr lt {p) O 



= Tr 



p rm 



(3.9) 



the semigroup j t is generated by (3.1), while T t by the dual equation: 

2n 



dOit) 
dt 



= i[H,0(t)}+ Yl C /H/ fatO(t)a v --{ata l/ ,0(t)}J . (3.10) 



Recalling the definition (2.9) and inserting in place of O suitable bilinears in the creation 
and annihilation operators, one sees that indeed the above equation readily implies (3.7). 

Being a first order differential equation involving finite-dimensional matrices, the evo- 
lution equation (3.7) can be easily integrated. Its solution involves the exponentiation of 
A and explicitly reads: 



V(£) = e* At • V(0) • e' A + 



tA 



f 

JO 



dr e TAt • B • e rA 



(3.11) 



In the next Sections, the dynamical evolution given by (3.11) will be used to investigate 
whether and under what conditions it can give rise to entanglement enhancement. 



4. ENVIRONMENT INDUCED ENTANGLEMENT GENERATION 

In order to study the possibility of entanglement generation by the external bath, it 
sufficies to consider a system formed by two oscillators. Henceforth, we specialize n = 2, 
so that the covariance V, as well as the Kossakowski C and the coefficients A, B, all 
become 4x4 matrices, and correspondingly r), a and A result 2x2 matrices. We shall 
concentrate our attention on discussing bath assisted entanglement production by purely 
dissipative mechanism; we shall therefore ignore any hamiltonian coupling between the 
two oscillators, setting in particular Coij = O.t 

t The entanglement power of purely hamiltonian couplings have been extensively stud- 
ied in the literature, e.g. see [6-10]. 
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Let us first analyze the possibility of entanglement creation at the beginning of the 
evolution, in the neighborhood of t = 0. As initial state of the two-oscillator system we 
shall choose a kernel p(z, z; 0) representing a separable Gaussian state and use the partial 
transposition criterion [42, 43] to check whether p(z,z;t) becomes entangled at a later 
time t. 

As well known, in quantum mechanics the operation of full transposition on any 
state p corresponds to the time-reversal transformation: it can be easily implemented 
in the holomorphic representation through the variable exchange Zi <-> z% in the kernel 
p(z,z). Recalling the expression in (2.4), this exchange equivalently corresponds to the 
transformations a <-> a* and f3 <-> f3 T on the submatrices defining the covariance G, and 
thus to the exchange ai <-> a\ in the definition (2.9) of symmetric covariance V. Clearly, 
the transformed V still satisfies the positivity condition (2.12) if the original covariance 
does, so that full transposition maps the set of Gaussian density matrices into itself. 

This is not the case for the operation of partial transposition, i.e. of transposition 
involving only one of the two oscillator system, say the first one, so that z\ <-> z\ . It clearly 
maps states into states for separable ones, but not in general for correlated ones: it thus 
provides a sufficient criterion for bipartite entanglement in any dimensions. 

For the case at hand, one finds the the operation of partial transposition with respect 
to the first system results in the following transformation of the 4x4 symmetric covariance: 



V-> V = T-V-T , T = 








1 







1 








1 











Vo 








1/ 



(4.1) 



We have seen that in order for a Gaussian kernel p(z, z) in (2.4) to represent a state, the 
condition (2.12) on its corresponding symmetric covariance V needs to be satisfied; if one 
further finds: 

V + |<0, (4.2) 

then the state is surely entangled. In the two-mode case we are studying, also the converse 
is true, namely, if V represents the symmetric covariance of an entangled Gaussian state, 
then (4.2) is necessarily satisfied [44]. 

In the case of the quasi-free markovian dynamics in (3.11), instead of dealing directly 
with the behaviour in time of V(£) and the inequality (4.2), we find it more convenient to 
consider the scalar quantity: 

C(f) = (*|(v(f) + !)|*) , (4.3) 

where is a four-dimensional complex vector, choosen in the null eigenspace of the 
matrix V(0) + S/2, so that 2(0) = 0. Then, the two oscillators, initially prepared in a 
separable Gaussian state, will start to become correlated by the noisy dynamics induced 
by the bath in which they are immersed if a suitable vector exists such that: 

d t Q(0)<0. (4.4) 
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By applying the partial transposition operation T to both sides of (3.7) and defining: 



|*) = T|#), E = T-E-T=( f" 3 ° ) , (4.5) 



a 3 

with o"3 the third Pauli matrix, the condition (4.4) can be cast in the following form: 

2<$|B|$) < <*|[A f •£ + £• A] |*) , (4.6) 

explicitly showing the dependence on the bath through the presence of the coefficients 
matrices A and B. 

As initial state for the two oscillators we take a separable state that is also pure; 
this is by no means a restriction: if the bath is not able to entangle pure states, it will 
surely not correlate their mixtures. Gaussian pure states are represented by holomorphic 
functions of the form (2.13), where now the matrix O is two-dimensional. Separability 
further impose the vanishing of the off-diagonal elements, so that O = diag(fii,fi 2 ), with 

< 1 to assure a finite state norm. Similarly, also the 2x2 submatrix a and /3, defining 
the covariances G and V, turn out to be diagonal: 



«=l iH o ir I , ■>=[ , ) • 



12 



o \ . o 



The condition (2.12), 



that assures the positivity of the corresponding Gaussian state, gives again the constraints 
< 1, as it should, since (4.8) is equivalent to norm convergence. 

The null subspace of the combination V + S/2 in (4.8) results two-dimensional; it is 
spanned by the vector |\&) of components (afl*, 60 2 , a, 6), with a and b abitrary complex 
parameters. Then, entanglement between the two oscillators will surely occur in baths for 
which (4.6) is satisfied with this choice of 

To show that indeed this is possible, let us initially prepare the two oscillators in their 
corresponding zero-temperature Fock vacua, so that fii = fi 2 = 0. By further choosing 
the parameters a and b real and equal, one finds that the inequality (4.6) reduces to: 

on + 022 < ^e(Ai2 + A 2 i) , (4.9) 

involving the entries of the 2x2 matrices a and A that appear in the Kossakowski matrix 
C in (3.6). As discussed in Section 3, the matrix C parametrizes the physical properties 
of the bath: it needs to be positive in order to comply with the requirement of complete 
positivity of the reduced dynamics of the two oscillators. Therefore, in order to be sure that 
baths satisfying (4.9) can actually be constructed, one needs to check the compatibility of 
(4.9) with the condition of positivity of the Kossakowski matrix. To explicitly show this, 
it is enough to take the matrices 17 and a diagonal and A real, with only the entry A21 non 
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vanishing. Then, the condition C > simply gives: fjn > 0, 1722 > 0, &n > 0, 022 > and 
A|i < ??22 Combining these conditions with that in (4.9), one obtains (for &u 7^ 0): 

l + ^<^< f?^V /2 . (4.10) 



en V a ii 

As a consequence, in order to be able to correlate the two oscillators, originally prepared 
in the separated Fock vacuum, it is sufficient to immerse them in a bath for which A21 is 
as in (4.10) and 7)22 greater than (an + d-22) 2 /&H- 

Environments that are able to entangle two oscillators that are prepared in separable 
temperature states can similarly be found. They are characterized by a Kossakowski matrix 
of the form: 

: rwi ; 



(4.11) 



where the parameters 77 and a are real and positive, while A is complex, satisfying the 
positivity condition: |A| 2 < r?cr. The initial state covariance is characterized by submatrices 
a and (3 as in (4.7); for simplicity, let us assume the single oscillator states to be equal, so 
that fii = O2 = ft- By taking the norm of the parameters a and b appearing in the vector 
|^) to satisfy \a\ = \b\ = l/y/2, and further adjusting their phases and that of A according 
to the relation Arg(A) = Arg(a) — Arg(6) = rr — Arg(O), one checks that the condition 
(4.6) assuring entanglement generation can be fulfilled by choosing 



M > -^17^ 1 (4-12) 



a — 


n\ 


V 


1 - 







which is compatible with the positivity condition of the above Kossakowski matrix provided 
we further take 77 > a. 



5. ASYMPTOTIC ENTANGLEMENT 



In the previous Section we have discussed examples of baths capable of correlate two 
independent oscillators through a purely noisy mechanism: no direct hamiltonian couplings 
between the two subsystems were present. This happens at the beginning of the evolution, 
as soon as t > 0. The test (4.4), on which this conclusion is based, is however unable to 
determine the fate of this entanglement as time becomes large. 

On general grounds, one expects that the effects of decoherence, counteracting entan- 
glement production, be dominant in the large time region, so that no quantum correlation 
is expected to be left at infinity. Nevertheless, there are situations in which entanglement 
is found to keep growing, reaching at the end an asymptotic non-vanishing value. An 
example is provided by baths for which the corresponding Kossakowski matrix takes the 
form in (4.11). 
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In this case, the evolution equation (3.1) can be rewritten in a simplified form by 
introducing the following set of independent oscillator variables: 

= oi + a 2 At _ a\ + 4 

^ f , (5.1) 

so that: [A, A^] = [B, B^] = 1 and all other commutators vanish. Then, taking into 
account (4.11), one finds that the equation (3.1) only involves the operators A and A^: 

^ = -io J [A^A,p\+L[p(t)} , (5.2) 

where, 

L\p]= V ([Ap,Ai] + [A,pA^)+a^p,A] + [A\ pA]) 

-\[A,[A,p]}-\*[A\[A\p}} . 

For sake of generality, we have also included an hamiltonian contribution; it comes from 
an effective Hamiltonian of the form (3.3), where, in analogy to (4.11), the 2x2 matrix Q 
is taken to be: 

u>ij = u J ^ , u > . (5.4) 

The general solution of the now single-oscillator dissipative evolution (5.2), (5.3), can 
be espressed in Gaussian form, using the holomorphic representation involving just one 
complex variable: A^ — > z, A — > d/dz. It can be explicitly expressed in terms of the 
analogs of the bilinears introduced in (2.6), now just numbers, a = (A 2 ) and (3 = (AA^), 
that parametrize the covariance of the corresponding Gaussian kernel solution p(z, z;t). 
Notice that we have (A) = (A^) = as a consequence of previous choice (aj) = (aj) = 0. 
One finds that the general solution of the evolution equations for a(t) and (3{t) derived 
from (5.2) is given by [45]: 

a(t) = e -2(.v-*+i*>)t ( ao _ ttoo) +aooi 
/ g(t)=e- 2 ^*(/3o-/3oc)+/3c 



'OO 1 



where the parameters «o = «(0) and (3q = (3(0) identify the initial values, while 

\*(<r-ri + iu) t] 
(rj — o) £ + uj z t] — a 

An asymptotic equilibrium state exists only when rj > a, condition that it will henceforth 
assumed. Then, the corresponding asymptotic Gaussian kernel is characterized by a co- 
variance with parameters a and (3 as in (5.6); it is not a thermal state, unless A = 0, in 
which case the temperature T is given by exp(u;/T) = rj/a.^ 



t Note that this result does not contradict the previous conclusions concerning bath- 
induced entanglement creation, that in fact requires A 7^ 0. 
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On the other hand, all averages involving B and result time independent, since 
(5.2) does not contain any dependence from these operators; for simplicity, we shall take 
them to be all zero, except for the bilinear (B B'): for reasons that it will appear clear 
soon, it is assumed to be equal to /?oo, the long-time asymptotic value of (AA^). 

Applying back the linear transformations (5.1), these results can now be used to study 
the evolution of the true physical oscillators, those described by canonical operators ai and 
a\. As a result, the two-oscillator state is represented by a time dependent Gaussian kernel 
p(z, z;t) giving rise to a covariance matrix V(t) as in (2.10): 

_ ( p(t) &(t) \ i 



Its time dependence is encoded into the two functions defined in (5.5), that appear in the 
entries of its two 2x2 submatrices etij and fay, explicitly, one finds: 



a(t) ( l l \ h ^ lfPW + Poc P(t)-Poo' 

2 




= -Yl Putt) = o I I • (5-8) 



P(t)~Poc P(t)+Pc 

For sake of simplicity, we shall further take ct(0) = «o = 0, so that V(0) only involves (3q 
and Pqq. The parameter (3 is not completely arbitrary: the positivity condition (2.12), 
V(0) + S/2 > 0, readily implies f3 > 1. Note that this two-mode initial state is mixed 
(compare it with the pure case in (4.7)) and further separable; indeed, after applying the 
partial transposition operation (4.1), one finds that 

V(0) + | > ; (5.9) 

indeed, this condition is equivalent to the inequality 2(3qPoo > Po + Poo, always satisfied 
for A), Poo > 1. 

We have previously shown that an environment described by a Kossakowski matrix in 
the form (4.11) is able to initially entangle two independent oscillators immersed in it. To 
see whether this still holds in the asymptotic long-time regime, one needs to examine the 
properties of the covariance of the corresponding equilibrium state Poo(z, z), obtained 
by letting a(t) — > a^, (3{t) — > 13^ in the expressions (5.7), (5.8). By using the partial 
trace criterion, will be entangled if and only if the 4x4 matrix Voo + S/2 posseses 
negative eigenvalues. In the present case, these eigenvalues can be analytically evaluated 
to be Poo ± y/AZ/2 - 1/2, A± = 1 + 2|a OG | 2 ± |aoo|(l + |ctoo| 2 ) 1/2 , for all four possible 
combinations of the ± signs. Recalling the definitions (5.6), one easily sees that the lowest 
eigenvalue — a/A + /2 — 1/2 can always be made negative by a suitable choice of the 
bath parameter A.t 

This result is remarkable: it shows that a dissipative quasi-free dynamics, generated 
by an equation of the form (3.1), (3.4), can produce quantum correlations even for large 



t The actual condition reads: |A| 2 > Arfa 2 [(i]-a) 2 +uj 2 ] / (if -a 2 ) 2 , which is compatible 
with the requirement of complete positivity, |A| 2 < rja, for sufficiently small ui. 
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times, allowing at the end an entangled equilibrium state. In view of the increasing interest 
that the theory of quantum information with continous variables is presently attracting, 
these results may be relevant both in phenomenological and experimental applications. In 
particular, the posssibility of mantaining bipartite entanglement in a noisy environment 
even for asymptotic long times may help the actual realizations of simple quantum devices 
in quantum optics and condensed matter physics. 
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